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ABSTRACT 

The photometric and spectroscopic observations of nearby type IIP supernova 
2004dj are presented. The ^^Ni mass estimated from the light curve is ~ 0.02 Mq. 
This estimate is found to be consistent with the Ha luminosity. SN2004dj reveals 
a strong asymmetry of the Ha emission hne at the nebular epoch with the shift 
of the maximum of —1600 km s~^. A similar asymmetric component is detected 
in H/3, [O I] and [Call] lines. The line asymmetry is interpreted as a result of 
the asymmetry of ^^Ni ejecta. The Ha profile and its evolution are reproduced 
in the model of the asymmetric bipolar ^^Ni and spherical hydrogen distribu- 
tions. The mass of the front ^^Ni jet is comparable to the central component 
and twice as larger compared to the rear ^^Ni jet. We note that the asymmetric 
bipolar structure of of ^^Ni ejecta is revealed also by SN1999em, another type 
IIP supernova. 



1. Introduction 



Supernova SN2004dj is discovered on 2004 July 31 in Sc galaxy NGC2403 (Nakano et 
al. 2004). According to spectra and light curve (Korcakova 2005) this is a usual type IIP 
supernova (i.e., with plateau in the light curve) detected in 30-45 days after the explosion. 
With the distance of 3.13 Mpc (Freedman et al. 2001) this is closest SNIIP after SN1987A. 
Given the peculiar nature of the latter SN2004dj turns out to be closest ever normal SNIIP 
observed. 

Type IIP supernovae are related with the explosion of a massive (10 — 25 Mq) red 
supergiant following the collapse of iron core. Three versions of the explosion mechanism 
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have been proposed: (1) neutrino mechanism (Colgate & White 1966; Buras et al. 2003); 
(2) magneto-rotational mechanism (Ardelyan et al. 2005); (3) rotational fragmentation of 
a new-born neutron star with the subsequent explosion of low-mass (~ 0.1 Mq) fragment 
(Blinnikov et al. 1990). Whatever happens the subsequent shock wave ejects the envelope 
with the sufficient internal energy to provide the SNIIP luminosity at the plateau stage 
(Grasberg et al. 1971). The shock wave propagation in the silicon mantle results in the 
synthesis of radioactive ^^Ni; its decay via ^^Ni-^^Co-^^Fe maintains the SNIIP luminosity 
at the tail of the light curve (Weaver & Woosley 1980). Neither of proposed explosion 
mechanism is able as yet to reproduce empirical values of the explosion energy and ejected 
s^Ni mass of SNIIP. 

The proximity of SN2005dj provides us an opportunity to study in greater details the 
observational display of an explosion of normal SNIIP. Unfortunately, supernova was dis- 
covered long after the outburst so the explosion energy cannot be recovered from a hydro 
modeling. Yet the light curve and Ha luminosity at the nebular epoch can provide us an 
estimate of the ^^Ni mass, another vital property of the explosion mechanism. 

Of particular interest would be the emission line profiles at the nebular epoch, which 
could provide us valuable information about the possible explosion asymmetry in SNIIP. 
The asymmetry issue emerged after the detection of the emission lines redshift in SN1987A 
(Phillips & Williams 1991) followed by the interpretation of this redshift as an effect of '"^^Ni 
ejecta asymmetry (Chugai 1991). A similar phenomenon was recovered recently in the type 
IIP supernova SN1999em (Elmhamdi et al. 2003a). Herant et al. (1992) conjecture that the 
^^Ni asymmetry as well as high pulsar velocities may be related to the violation of a point 
symmetry in core-collapse supernovae owing to the large scale neutrino convection. 

Below we present results of spectral and photometric observations of SN2004dj in Special 
astrophysical observatory (SAO) at Zelenchuk. Apart from a qualitative analysis of the 
observational data we present estimates of ^^Ni mass in SN2004dj from the light curve and 
the Ha luminosity. The observed spectra reveal unusually strong asymmetry of the Ha 
emission line. Given the importance of this phenomenon we concentrate on the study of the 
line asymmetry and its relationship with the asymmetry of ^^Ni ejecta. 

Henceforth we use the distance D = 3.13 Mpc (distance module // = 27.48) according to 
cepheids (Preedman et al. 2001) and the redshift of 129 km s~^ according to LEDA database. 
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2. Observations 

Broad band BVRc photometry of SN2004dj was carried out at 1-m Zeiss telescope 
of SAO using CCD Electronika K-585 and EEV 42-40. The photometry data were also 
obtained at 60-cm telescope of Crimean station of SAI with CCD VersArray. These three 
detectors are shortened as K, E, and VA respectively. A local photometric standard around 
supernova is based on the U BVR standard in the vicinity of NGC 2403 (Zickgraf et al. 1990) . 
The reduction of CCD images is made with the WinFITS package (V.P. Goranskij) using a 
method of corrected aperture measurements that permits us to deal with the inhomogeneous 
background in the vicinity of supernova. The accuracy of the photometric data is 0.005-0.01 
mag. The data are presented in Table 1. The columns in order are Julian dates, magnitudes 
and the detector type. The first line of the Table 1 contains the photometry data on 2001 
January 19, long before the explosion. We believe that these magnitudes refer to the star 
cluster S96 on which SN2004dj is superimposed. 

The log of spectral observations which includes observation moments, spectrograph type, 
spectral range and resolution is presented in Table 2. Eight spectra were taken between 2004 
October 18 and 2005 June 9. The first spectrum is obtained at the transition between the 
plateau and tail, whereas the rest ones were taken at the light curve tail. Three types of 
spectrograph were used in observations: (i) Scorpio, the universal focal reducer (Afanasiev 
& Moiseev 2005) of 6-m large azimuthal telescope (LAT), multi-pupil fiber spectrograph 
(MPFS) of LAT, and UAGS spectrograph of 1-m telescope of SAO. CCD images of spectra 
taken at Scorpio and UAGS are reduced using MIDAS package with the usual procedures 
of the extraction of cosmic rays, field fiattening, and spectra extraction. The reduction of 
spectra taken at MPFS is performed using package written by V.L. Afanasiev, A.V. Moiseev 
and P.K. Abolmasov. In two cases (spectra No. 2 and 6 in Table 2) the atmospheric 
dispersion strongly modifies the spectral energy distribution. These spectra will not be used 
for the Ha flux measurements, although they may well be used for the line proflle analysis. 



3. Photometry and Ni mass 

The V light curve according to our data and other available data^ is plotted in Fig. 1. 
The time is counted from the adopted explosion moment JD2453170 (i.e., 2004 June 13). This 
choice is based upon the assumption that light curves of SN2004dj and SN1999gi, standard 
SNIIP (Leonard et al. 2002), are similar. This assumption is supported by the fact that the 



^ http: / /www. astrosurf . com / snweb2 /2004dj /Meas .htm 
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same fitting curve describes both sets of data reasonably well (Fig. 1). The approximation 
suggests the exponential behavior of the flux in V band at the tail F{V) oc exp(— i/111.26d) 
in accordance with the radioactive decay of ^^Co-^^Fe. 

The two last measurements (day 266 and 267) show excess ^ 0.17 mag over the expo- 
nential behavior (Fig. 1). Most of this excess is hkely related to the contribution of the 

star cluster S96. Given the prc-cxplosion magnitude (Table 1) we find that the star cluster 
on day 266.5 must contribute 0.14 mag assuming exponential behavior of the intrinsic V 
magnitude of SN2004dj. The light of S96, indeed, accounts for, practically, all the excess. 

Galactic reddening towards NGG2403 is E{B - V) ^ 0.062 (Schlegel et al. 1998). 
This value is adopted as a total reddening of supernova in our and host galaxy. With this 
reddening the intrinsic B — V color of SN2004dj according to our and other available data (see 
previous footnote) is displayed in Fig. 2. We plotted also the SN1987A color (Catchpole 
et al. 1987, 1988) and the lower envelope of intrinsic colors of SNIIP (Elmhamdi et al. 
2003b). The B - V color of SN2004dj coincides wih the lower envelope of SNIIP colors, 
which indicates that the adopted reddening is close to the actual value. On the other hand 
Wang et al. (2005) upon the basis of the population synthesis of the star cluster S96 derive 
the reddening E{B — y) ~ 0.35. The contradiction may be resolved if the light form S96 is 
absorbed in the cool shell around this cluster, while the supernova resides outside the shell, 
closer to the observer. 

The V - R color of SN2004dj (Fig. 2, inset) around day 150 is by 0.3 mag lower than 
that of SN1987A. However, the color difference gets smaller with time and disappears around 
day 200. This indicates that the spectral energy distributions (SED) for both supernova in 
the long wavelength range (A > 5000 A) at t > 200 d are similar. This fact is of importance 
for the photometric estimate of ^^Ni mass in SN2004dj. 

The amount of ^^Ni in the envelope can be found from a comparison of My absolute 
magnitudes of SN2004dj and SN1987A at the early nebular epoch (t < 250 d) when the enve- 
lope is optically thick for gamma-rays, so one may adopt that all the radioactive luminosity 
is instantly re-emitted in the optical. The method assumes that the distance and reddening 
are well known and the bolometric correction is similar for both supernovae. The latter 
seems to disagree with the apparent difference of their B — V colors. However, one should 
keep in mind that at this epoch the bulk of the energy is emitted in the long wavelength 
range (A > 5000 A), in which SED are similar around and after day 200 as indicated by 
the their V — R colors (Fig. 2). Therefore, using SN1987A as a template we hopefully can 
recover the amount of ^^N in SN2004dj from V magnitude around day 200. 

Using the exponential approximation of the V light curve of SN2004dj (Fig. 1) we get 
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V = 15.32 on day 200. With the reddening E{B -V) ^ 0.062 and standard extinction 
law (Seaton 1979) the absorption is Ay = 0.19. With the distance module fi = 27.48 for 
NGC2403 one gets My = -14.01 on day 200. On the other hand, for SN1987A at this 
epoch V = 5.09 (Catchpole et al. 1988). Adopting the distance 50 kpc and absorption 
Av = 0.6 one gets My = -14.01 on day 200. Thus on day 200 SN1987A is brighter in V 
band compared to SN2004dj by 1.53 mag. With the ^^Ni mass 0.075 M© in the SN1987A 
envelope (Suntzeff & Bouchet 1991) we conclude that the amount of ^^Ni in SN2004dj is 
M(56Ni) ^ 0.02 Mq. 

The uncertainty of the "''^Ni mass related to the error in distance (Freedman et al. 
2001) is of 20%. The uncertainty related with the extinction unlikely exceeds 10% taking 
into account low reddening indicated by the extremely "blue" color B — V oi SN2004dj. 
Some contribution to the error of ^^Ni mass may be related to the possible difference of 
bolometric corrections of SN2004dj and SN1987A. To summarize, the total error of the ^^Ni 
mass estimate is ~ 30%. 

Independently, the ^^Ni mass may be estimated using the empirical relation between 
amount of ^^Ni and parameter S — dV/dt that characterizes the steepness of the V light 
curve at the transition between plateau and tail (Elmhamdi et al. 2003b). In the case of 
SN2004dj we find S = 0.17±0.04 mag. day~^. Using the empirical relation from (Elmhamdi 
et al. 2003b) one gets M(^^Ni) = 0.013 ± 0.004 Mq that is consistent within errors with the 
estimate found above. 



4. Spectrum and its evolution 

4.1. Qualitative analysis 

The observed spectra of SN2004dj are presented in Fig. 3. The last two spectra (Table 
2) are merged into the single spectrum since at this epoch the spectral evolution is slow. The 
appearance and evolution of SN2004dj spectrum are similar to those of ordinary SNIIP at 
the nebular epoch. The late time spectrum is dominated by emission lines of Ha, [O I] 6300, 
6364 A and [Call] 7291, 7324 A on the background of the relatively weak quasicontinuum. 
According to general wisdom, the latter is composed by the emission of numerous metal lines, 
mostly Fell, that originates as a result of the absorption of ultraviolet photons generated 
in central parts of the envelope due to the reprocessing of the radioactive luminosity (Xu & 
McCray 1991). The scattering and absorption of the quasicontinuum in metal and hydrogen 
lines in outer layers produces broad P Cygni profiles, e.g., H/3, Fell lines of the multiplet 
42, Nal doublet (Fig. 3). 
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A spectacular feature of the SN2004dj spectra, never observed before in SNIIP, is the 
strong blueshift of emission hne maxima, most apparent in Ha (Fig. 4). The profile of this 
line is clearly asymmetric with the maximum at —1600 km s~^ in early (127 and 131 days) 
spectra. The blueshift decreases with time down to -600 km s^^ on day 330. The blue peak, 
although weaker, is present also in B.f3. Remarkably, on day 131 the Ha profile shows double 
horn structure with the red horn shift of ~ 1600 km s~^. This indicates the bipolar, although 
asymmetric, structure of the hydrogen excitation. 

Among type HP supernovae only two cases of significant Ha line asymmetry were men- 
tioned earher: SN1987A (Phillips & Williams 1991) and SN1999em (Elmhamdi et al. 2003a). 
In both cases the redshift was reported. The redshift in SN1987A has been interpreted as 
a result of the ejection of significant amount of ^^Ni in the far hemisphere in otherwise 
spherical envelope (Chugai 1991). The similar interpretation was proposed for the redshift 
of Ha and Hel 10830 A fines in SN1999em (Elmhamdi et al. 2003a). We adopt the same 
explanation for the Ha asymmetry in SN2004dj with the only difference that in this case the 
^^Ni distribution is skewed towards the observer. 

The manifestation of the fine asymmetry in SNHP is displayed in Fig. 5, which shows Ha 
and [CaH] doublet in SN2004dj on day 149, SN1987A on day 198 (Phillips & Wilfiams 1991) 

and SN1999em on day 140 (Elmhamdi et al. 2003a). Apart from the sign of asymmetry, 
SN2004dj differ by the larger shift and the amplitude of the asymmetric Ha component. The 
[Call] doublet of SN2004dj differs also very much from that in other supernovae (Fig. 5). 
SN2004dj shows the striking blue peak in [Call] at —1600 km s~^ that is seen in Ha. This 
peak reflects, presumably, the local overexcitation of Ca II related to the asymmetry of ^^Ni 
ejecta likewise in Ha. 

Other spectral differences between SN2004dj and SN1999em related to different signs 
of asymmetry are noteworthy. The absorption component of Ha in SN2004dj is markedly 
weaker than in SN1999em (Fig. 6). This is caused by the presence of the blue emission 
peak in SN2004dj that fills in the absorption. Absorption component of H/3, on the contrary, 
is stronger in SN2004dj. The reason for that is also the ejection of ^^Ni primarily towards 
the observer. This causes higher hydrogen excitation in the near hemisphere of SN2004dj 
compared to SN1999em where most of ^^Ni is shifted towards the far hemisphere. As a result 
the H/3 optical depth in the near hemisphere of SN2004dj is greater than in SN1999em. 

Interestingly, spectra of SN2004dj at the plateau epoch (Korcakova et al. 2005) demon- 
strate larger expansion velocities of the Ha absorption minimum compared to SN1999em. 
On JD2453221 (day 51 for the assumed outburst date) radial velocity of the Ha absorption 
minimum is -6955 km s~^ (Korcakova et al. 2005) compared to -5000 km s~^ in SN1999em 
on 54 day (Elmhamdi et al. 2003a). We believe that higher velocity of SN2004dj refiects 
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the higher degree of the hydrogen excitation in the near hemisphere rather then higher 
energy-to-mass ratio. 

Another remarkable pecuharity of SN2004dj is the quasicontinuum shape in the range 
of 5000-5700 A which differ from that in SN1999em. Note, the quasicontinua of SN1999em 
and SN1987A are rather similar in the same range. The plausible explanation involves the 
different sign of ^^Ni asymmetry in SN2004dj and SN1999em. Actually, the ^^Ni displace- 
ment should result in the different conditions for the quasicontinua formation in asymmetric 
hemispheres because of the different excitation degree. 

4.2. Ha luminosity and ^^Ni meiss 

The integrated Ha flux is measured in the spectra which were cahbrated using V and 

Rc photometric data. The resulting Ha flux is the average between fluxes derived from V 
and Rc calibrations. We omit spectra on days 131 and 237 which are strongly affected by the 
atmospheric dispersion, and the spectrum on day 330 which cannot be calibrated confidently 
because the strong [Call] emission is missing. The obtained Ha fiuxes, uncorrected (F) and 
corrected (^F'^°") for the dust absorption, are given in Table 3. The last column is the Ha 
luminosity. The error of the flux (0.12 dex) is estimated as half of the maximal difference of 
ffuxes derived from V and Rc calibrations. 

The Ha luminosity at the nebular epoch is used here to estimate ^^Ni in SN2004dj 
likewise this was done earher for SNHP (Elmhamdi et al. 2003b). The method suggests the 
computation of Ha luminosity powered by the radioactive decay ^^Co - ^^Fe. Three levels 
plus continuum hydrogen atom is adopted, although recombinations on the excited levels 
are included. The photoionization from the second and third level by continua is taken into 
account. Radiation transfer in continua is treated in the escape probability approximation 
(Chugai 1987). The two-zone model envelope consists of the macroscopically mixed core and 
outer envelope. The model parameters are the mass M, kinetic energy E, ^^Ni mass, mass 
fraction in the core /mix and the fraction of the core residing in metals and helium f^. The 
quasicontinuum SED is described by the plankian spectrum with the color temperature T^. 
In contrast to the previous version here the electron temperature is calculated from the 
energy balance which takes into account cooling in CaH, MgH, FeH (approximately), 01 
lines, //-emission and adiabatic cooling. We abandon the approximation of the homogeneous 
density assuming that the density in the core is three times larger than the outer shell density. 
This mimics the density distribution in the real SNHP envelope. 

The model takes into account the outcome of the macroscopic mixing of ^^Ni in the 
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exploding star with the stratified chemical composition. We adopt that fragments of the 
^^Ni material are embedded in cocoons composed of metals and helium, i.e., of the material 
lacking the hydrogen. The optical depth of the cocoon for gamma-rays is determined by 
the mass and the surface area of cocoons, = ^/m/mix^/ (47r/s-R^), where k = 0.03(1 + X) 
cm^ g-i (Fransson and Chevalier 1989), X is the hydrogen abundance, Ri is the radius of 
the mixed core, fs is the mixing parameter. The latter is the ratio of the cocoon surface 
area to AnRl. Below we show results for fs — 2 (moderate mixing) and /g = 20 (strong 
mixing). The adopted ejecta mass is M = 13 Mq, presumably the mass of SN1999em 
ejecta (Nadyozhin 2003) that corresponds to the 15 Mq main sequence mass. The kinetic 
energy is = 10^^ erg, the typical value for SNIIP (Nadyozhin 2003). Other parameters 
are /mix = 0.4, /m = 0.3, i.e., ~ 1.56 A/q of the core resides in metals and helium. The 
remainder of helium, from several tens to ~ 1 Mq is presumably mixed with the hydrogen. 
Note, for the star with the initial mass M «i 15 Mq the predicted mass of newly synthesised 
metals and helium in the supernova is about 3 Mq (Woosley & Weaver 1995), in agreement 
with our assumptions. 

The optimal models of the Ha luminosity in SN1999em at the nebular epoch t > 130 
days are presented in Fig. 7 for the two values of the mixing parameter /g. The adopted 
^^Ni mass M(^^Ni) = 0.027 Mq is equal to the value recovered from the tail of V light curve 
(Elmhamdi et al. 2003b). The color temperature of the quasicontinuum is the only free 
parameter. Its value turns out to be — 5500 K for the model with /s = 2 and — 5000 K 
for /, = 20. 

For SN2004dj we adopt the same model as for SN1999em. The only difference is the 
^^Ni mass, which is M(^^Ni) = 0.02 Mq, the value recovered from the V light curve. The 
plot (Fig. 7, lower panel) shows that the model of the Ha luminosity with this amount of 
^^Ni fits to the observational data within errors. Of course, this fact cannot be considered 
as a proof that the ^^Ni mass is actually ~ 0.02 Mq given uncertainties of other parameters. 
On the other hand, this indicates that the ^^Ni mass recovered from the V light curve is 
approximately consistent with the Ha luminosity. 



5. Line asymmetry and ^^Ni distribution 

The line profile formed in the supernova envelope with the kinematic of the homologous 
expansion {v = r/t) depends on the distribution of the level populations and the continuum 
optical depth. At the nebular epoch the absorption in the quasicontinuum is weak in the red 
part of the spectrum (A > 5000 A). The profile of emission lines of Ha, [01] 6300, 6364 A 
and [Call] 7291, 7324 A thus contain the information about the asymmetry of the radiation 
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sources throughout the envelope including deep central zone. 

5.1. [OI] 6300, 6364 A lines 

The lines in [01] and [Call] doublets are partially superimposed, so it would be rea- 
sonable to recover true line profile using representation of each line as a superposition of the 
minimal number of gaussian components. Each gaussian component is characterized by its 
radial velocity v^, Doppler width w and relative amplitude A. 

Oxygen doublet [O I] 6300, 6364 A is a convenient case for the decomposition because the 
velocity interval between lines is large (3000 km s~^). Using only two gaussian components, 
central and blueshifted, we are able to reproduce successfully the doublet profile on day 237 
and 330 (Fig. 8). Parameters of components (displacement, width, amplitude) are given in 
Table 4). The line profiles of [O I] and Ha on day 237 are quite similar (Fig. 4); on day 330, 
however, the difference is remarkable: the [O I] maximum is at zero, while the Ha maximum 
is blueshifted at ^ —700 km s~^. This fact reflects, probably, the different distribution of 
hydrogen and oxygen. 

The increasing role of the central component can be explained in a picture of the two 
^^Ni components. Let the central component be embedded into the cocoon that is optically 
thick for gamma-rays and does not contribute in the [O I] doublet emission. Note, the cocoon 
material may be even composed by oxygen and yet be unseen in [O I] lines. This is the case, 
if formed CO and/or SiO molecules cools the oxygen material (Liu & Dalgarno 1995). As 
the envelope expands, the cocoon gets optically thin and gamma-rays efficiently excite outer 
oxygen material thus increasing the role of the central component. A similar effect in Ha 
will be studied below (Section 5.3). 

For the doublet decomposition we used some value of the doublet ratio R — /(6300) / 1(6364) . 
The latter is equal to the ratio of line specific luminosities hvikUkAkiPik [erg s~^ cm~^], where 
Pik is the escape probability, while the other notations are obvious. The optimal R values 
(Table 4, last column) are lower than nebular one, R = 3, which implies that [01] lines are 
optically thick. The similar situation takes place in other SNIIP (Spyromilio & Pinto 1991; 
Chugai 1988). The doublet ratio is determined primarily by the oxygen concentration and 
only weakly depends on the electron temperature, since the excitation potential of the lower 
level of the transition of 6364 A is only 0.02 eV. Since the oxygen is mostly neutral the ratio 
R permits us to recover the oxygen density in the line forming zone. The observed ratio 
R = 1.33 on day 237 leads to the oxygen number density n(0) = 6.63 x 10*^(237 d/t)^ cm-^ 
where t in days and — 5000 K is adopted. The density evolution n(0) oc predicts 
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the ratio R — 1.77 on day 330, in good agreement with the value R — 1.8 found from ob- 
servations (Table 4). This confirms the adopted explosion date of SN2004dj, although the 
accuracy of the reconstruction of the explosion date from the ratio R is not better than ±30 
days. 

5.2. [Call] 7391, 7324 A lines 

The lines of [Call] doublet also can be represented by two gaussian components, central 
and blucshifted (Fig. 9, Table 5). However, the broad central component shows some 
deviation from gaussian shape and requires an additional weak narrow component (Fig. 9, 
inset). On day 237 when [Call] and [01] doublets can be compared, the line profiles are 
roughly similar, although some differences arc apparent. Particularly, the blue component 
of [Call] shows the same shift but 1.6 times more narrow. The central component, on the 
contrary, is 1.7 times broader than in [01] case. The average radial velocity of [Call] line 
(e.g., 7291 A) is -150 km s~^ on day 237, i.e., lower than that of [01] line at the same epoch 
(-770 km s-i). 

These differences are likely caused by the formation of these lines in different zones. 
Actually, stellar evolution for massive stars (e.g., 15 Mq) and supernova models predicts 
that Ca is synthesised in two processes: static oxygen burning in the presupernova and 
explosive oxygen burning during the shock wave propagation (Woosley Sz Weaver 1995). 
Calcium, therefore, is expected to reside in products of the oxygen burning (in Si/S matter), 
not in the oxygen matter. Moreover, the [Call] 7391, 7324 A lines may also form in hydrogen 
or helium material (Li & McCray 1993) with the main sequence abundance. 

Given uncertainty of the site of the [Call] doublet formation, it is interesting that the 
doublet ratio on day 149 is R = J(7291)//(7324) ^ 1.1, i.e., lower than the nebular one 
{R = 1.5). If this is caused by the optical depth effect, then the required number density of 
Ca should be n(Ca) > 3.6 x 10^(149 d/t)"^ cm^"^. Assuming that the [Call] doublet is emitted 
by the hydrogen matter we find the required density of hydrogen for normal Ca abundance 
p > 3.7 X 10~^^(149 d/t)^ g cm~^. This should be compared with the oxygen density found 
from oxygen doublet ratio, p 7x 10~^^(149 d/t)^ g cm~^. The latter turns out to be five time 
lower compared to the required hydrogen density, quite reverse to expectations. Actually, 
the analysis of SN1987A shows that the oxygen is distributed in dense condensations with 
the filling factor of / ~ 0.1 and the oxygen density is by order of magnitude greater than the 
average density (Spyromilio & Pinto 1991; Chugai 1988). The found controversy indicates 
that the hypothesis of the [Call] doublet formation in the hydrogen material faces serious 
problems. However, given uncertainty of the quasicontinuum behavior in the 7300 A band 
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one cannot preclude that we underestimated the doublet ratio. We thus do not rule out 
the formation of [Call] doublet in the hydrogen material although this possibility seems 
doubtful. 



5.3. Modeling Ha asymmetry 

Here we address the Ha hne profile and its connection with the asymmetry of the ^^Ni 
distribution. We facilitate the physics of the energy cascade resulting in the Ha emission 

and assume that the specific Ha luminosity, 1] = hu23n^,A;i2l323, is proportional to the local 
deposition rate rj = Ce^ with the constant factor C. To obtain the Ha line profile one 
has to define the ^^Ni distribution and to calculate the distribution of the deposition e^. 
The gamma-ray transfer is computed in a single flight approximation with the absorption 
coefficient adopted above (section 4.2). 

The model envelope expands homologously {v = r/t) and has the spherical mass distri- 
bution (but ^^Ni). The density is described by the exponential distribution p = po cxp(— f /fo) 
with Po t~^. The exponential law is fairly good for the bulk of the inner envelope where the 
power index \dln p/dlnv\ gradually increases outwards (e.g., Utrobin 2004). The parameters 
Po and Vq are determined by M and E. We adopt as above M = 13 Mq and E = 10^^ erg. 

The ^^Ni distribution is set by three components: central homogeneous sphere and two 
cylindrical jets with the angle 6 between the jet axis and the line of sight (Fig. 10). In 
the regions, where jets are superimposed on the central spherical component, the density 
is assumed to be equal to the jet density. We found that the opposite condition is less 
consistent with observations. The density of ^^Ni in jets is assumed to decrease hnearly 
outwards down to zero at the outer jet boundary. The chemical composition of the envelope 
is homogeneous with the hydrogen abundance X = 0.7. In the central zone of the radius 
Rs = Vgt the matter distribution is presumably inhomogeneous with the filling factor of the 
hydrogen matter linearly decreasing inwards from unity at Rs to zero in the center. This 
assumption is consistent with intuitive picture of the mixing of helium and metals with the 
hydrogen envelope. 

The study of the parameter space of the model led us to conclude that only the version 
with cocoons around ^^Ni material is able to reproduce the evolution of Ha profile. The 
cocoon model has been already used above to describe the evolution of Ha luminosity. Here 
this model should be refined given more complicated geometry of mixing. We assume that 
^^Ni is represented by an ensemble of spherical fragments each coated by cocoon of the 
optical depth Tc for gamma-rays. The average gamma-ray emissivity is then multiplied by 
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the attenuation factor exp(— Tc). We adopt Tc = const in the central zone of the radius Rs, 
while in jets outside the central zone, the cocoon optical depth decreases with the velocity 



where is the outer velocity A;-th jet {k = 1, 2). The adopted power index is p — 3. The 
parameter tq may be written via the cocoon mass of as Tq — kMf./ AirfgRl where /g is the 
mixing parameter introduced earlier (Section 4.2). The decrease of Tc outwards corresponds 
to the increase of the mixing degree. This is just what is expected from the mixing physics: 
the material propagated through the larger mass should experience larger mixing degree. 

Parameters of the optimal model, i.e., the angle 9 between the hne of sight and the axis 
of ^^Ni, the expansion velocity of the central component Vs, boundary velocity of jets {vi, and 
V2), the mass of jets, and cocoon mass for the mixing parameter /s = 2 are given in Table 5. 
The accuracy of the determination of velocities is about 20%, angle 9 is determined with the 
accuracy of 10°. Interestingly, while the model of the Ha luminosity does not impose strict 
constraints on /g, we find that assuming a cocoon mass Mf. ~ 1 Mq requires a moderate 
mixing, /s ~ 2. 

Model profiles of Ha with and without cocoon are plotted along with observations in 
Fig. 11. The model with cocoon is apparently more adequately describes profile and its 
evolution during 131-330 days. Note, within the same model wc arc able to describe the 
non-trivial change from asymmetric two-horn profile on day 131, through the plateau plus 
blue peak on day 149, towards smooth asymmetric profile with the blueshifted maximum. 

The specific features of the recovered ^^Ni distribution are the bipolar structure and 
the deviation from the point symmetry. The mass of the front ^^Ni jet is comparable to the 
central component and twice as larger compared to the rear ^^Ni jet (Table 5). Note, there 
is no need to admit the asymmetry of hydrogen distribution, which otherwise would require 
extremely large degree of the supernova asymmetry. 

The one-dimension projection of the recovered ^^Ni distribution on the line of sight and 
the jet axis are shown in Fig. 12. For comparison we show also the recovered hne profile 
of [Call] 7291 A on day 149 (Fig. 9), when the gamma-ray transfer effect is small and, 
therefore, the sources in [Call] lines follow the distribution of ^^Ni. While there is a general 
correspondence between [Call] hne profile and ^^Ni distribution we see apparent differences. 
The ^^Ni projection on the line of sight does not show strong central peak that is present in 
[Call] line. This disparity indicates that the Ca distribution deviates from ^^Ni distribution. 

Interestingly, we found that the major mass of ^^Ni should be isolated from the hydro- 
gen by the cocoon that does not contain any hydrogen. Moreover, the mixing degree that 



V as 
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determine the the thickness of cocoon should not be large. The similar property of the ^^Ni 
distribution has been revealed by Liu & Dalgarno (1995) in SN 1987A. To account for the 
evolution of CO 2.3 mem emission they introduced a variable reduction factor for gamma- 
rays that decreases with time. This factor, we believe, reflects the existence of cocoons 
around ^^Ni fragments in SN1987A as well. 

6. Discussion and conclusion 

We have presented results of the photometry and spectroscopy of SN2004dj, the usual 
type IIP supernova. We have estimated ^^Ni mass from the comparison of light curves of 
SN2004dj and SN1987A and from the modeling of Ha luminosity. Both estimates suggests 
the ^^Ni mass of ~ 0.02 Mq. This is 3.8 times less than in SN1987A and comparable with 
s^Ni mass in SN1999em (0.02 - 0.027 M©, Elmhamdi et al. 2003a, 2003b). 

The principal result of spectroscopic study of SN2004dj is the detection of a strong 
blueshift of the Ha line that is interpreted as a result of the asymmetric ejection of ^^Ni in 
the spherically-symmetric envelope. The modeling of the Ha line profile led us to conclude 
that the structure of the profile and its evolution are well reproduced in the model of the 
asymmetric bipolar ^^Ni ejecta with the two-fold difference of jet masses and the comparable 
masses of the massive jet and the central component. 

We have found signatures of the bipolar structure in Ha and [Call] line of SN1999em, 
which indicate that the asymmetric bipolar ^^Ni ejecta is not exception for SNIIP. Although 
in SN1987A spectral signatures of the bipolar structure are absent, Wang et al. (2002) see 
bipolar structure of inner region of SN1987A in images taken by HST. The detection of 
bipolar structure of ^^Ni ejecta at least in two normal SNIIP implies that bipolar explosion 
of the core-collapse supernovae can occure for a moderately rotating red supergiant core. 

The significant deviation from the point symmetry in the ^^Ni distribution revealed by 
at least three supernovae, namely SN1987A, SN1999em and SN2004dj, suggests that the 
total momentum may be essentially non-zero for synthesised iron-peak ejecta in SNIIP. In 
this case the finite momentum is presumably compensated by the envelope or/and neutron 
star. According to recent study of proper motions of pulsars (Hobbs et al. 2005) neutron 
stars at birth acquire high velocities about 400 km s^^ in average. These velocities unlikely 
could be explained by the sling effect of supernova in close binary, so the kick is the preferred 
explanation. The asymmetry of ^^Ni ejecta in SNIIP combined with high pulsar velocities 
is thus a strong argument that the explosion of core-collapse supernovae is accompanied by 
the violation of the point symmetry. 
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At present the likely mechanism for the asymmetry in core-collapse supernovae is the 
large-scale neutrino convection (Herant et al. 1992; Scheck et al. 2004), although the ejection 
of the supernova envelope with the required kinetic energy of ~ 10^^ erg remains problematic 
(Buras et al. 2003). 

The detection of the bipolar structure of the ^^Ni ejecta in normal SNIIP raises another 
intriguing question, whether the bipolar structure could emerge in the neutrino-driven explo- 
sion mechanism? This problem was recently discussed in connection with the interpretation 
of the bipolar ejecta of Cas A and SN1987A (Burrows et al. 2005). Authors considered 
the role of the moderate rotation in the creation of jets within neutrino- driven mechanism. 
However, it should be emphasised, that both these supernovae might in principle originate 
from presupernovae with the relatively fast rotation: Cas A is SNIb/c, and the nebular ring 
around SN1987A indicates fast rotation too. As to SN2004dj and SN1999em, which are 
normal SNIIP, here we deal with cores of red supergiants that are presumably slow rotators. 
Therefore, it might be well that the origin of the asymmetric bipolar ^^Ni ejection in SNIIP 
has nothing to do with the rotation. 

We are grateful to V.L. Afanasiev, N.B. Borisov nad S.S. Kaisin for the help in ob- 
servations. The work is partially supported by RFBR grants 4-02-17255, 03-02-16341 and 
04-2-16349. 
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Table 1: SN2004dj photometric data 
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Table 2: Log of spectral observations 
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Table 3: Flux and luminosity of Ha 
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Table 4: Parameters of components of lines of [01] and [Call] doublets 
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Table 5: Parameters of ^^Ni distribution 
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Fig. 1. — Light curves in V band of SN2004dj and SN1999gi. Triangles are our data, points 
are amateurs data. Solid line shows the analytical fit, the same in both cases. The time is 
counted from the adopted explosion moment (2004 June 13). 
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Fig. 2. — Color B — V oi SN2004dj. Triangles are our data, points are amateurs data. 
Solid line is the SN1987A color, dashed line is the lower envelope oi B — V colors of SNIIP 
(Elmhamdi et al. 2003b). Inset shows the V — R color of SN2004dj (crosses) in comparison 
with the color of SN1987A (line). 
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Fig. 3. — Spectra of SN2004dj at different age and major identifications. The wavelenth is 
corrected for the galaxy redshift. The age next to spectrum is measured from the explosion 
(2004 June 13). 
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Fig. 4. — Evolution of Ha in spectra of SN2004dj. The age is indicated next to line of a 
corresponding width. 
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Fig. 5. — Ha line and [Call] doublet in spectrum of SN2004dj on day 149, in spectrum of 
SN1987A on day 198 and in spectrum of SN1999em on day 140. The radial velocity of [Call] 
is given for the line of 7291 A. 
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Fig. 6. — Spectra of SN2004dj and SN1999em at similar phases. Note, the difference of 
profiles and quasicontinua around 5600 A. 
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Fig. 7.— Evolution models of Ha for SN1999cm and SN2004dj. Data for SN1999em (upper 
panel, squares) are taken from Elmhamdi (2003b). Model of SN1999em suggests ^^Ni mass 
of 0.027 Mq and two options of mixing parameter /s = 2 {thick line) and fs — 20 {thin line). 
For SN2004dj the model is the same as for SN1999em, but ^^Ni mass is 0.02 Mq. 
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Fig. 8.— Reconstruction of the line profile of [O I] 6300, 6364 A doublet in SN2004dj on day 
237 and 330. Thick line is the model, thin line is observation. Inset shows model line profile 
{thick line) along with gaussian components. 
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Fig. 9.— Reconstruction of the line profile of [Call] 7291, 7324 A doublet in SN2004dj on 
day 149 and 237. Thick line is the model, thin line is observation. Inset shows model line 
profile {thick line) along with gaussian components. 
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Fig. 10. — Schematic representation of the spherical envelope of SN2004dj with bipolar Ni 
distribution, The ^^Ni distribution includes the central component S and two jets Jl and J2. 
Arrow shows the line of sight. 
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Fig. 11. — Models of the Ha in SN2004dj with the asymmetric bipolar ^^Ni distribution. 
Thick line is the model, thin line is observations. Models of ^^Ni distribution without cocoon 
{left panels) and with cocoon {right panels) are shown. 
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Fig. 12. — Mass distribution of ^^Ni along the radial velocity. Upper panel shows projection 
on the line of sight {thick line) together with the model profile of [Call] 7291 A line on day 
149 (Fig. 9). The lower panel shows the projection of ^^Ni distribution on the polar axis of 
ejecta. 



